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We study the sizes and thermal properties of glueballs in a three dimensional compact Abelian 
gauge model on improved lattice. We predict the radii of ~ 0.60 and ~ 1.12 in the units of string 
tension, or ~ 0.28 and ~ 0.52 fm, for the scalar and tensor glueballs, respectively. We perform a well 
controlled extrapolation of the radii to the continuum limit and observe that our results agree with 
the predicted values. Using Monte Carlo simulations, we extract the pole-mass of the lowest scalar 
and tensor glueballs from the temporal correlators at finite temperature. We see a clear evidence 
of the deconfined phase, and the transition appears to be similar to that of the two-dimensional 
XY model as expected from universality arguments. Our results show no significant changes in the 
glueball wave functions and masses in the deconfined phase. 
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I. INTRODUCTION 

The prediction of glueball masses has longbeen at- 
tempted in lattice gauge theory calculations [HM & 0> • 
These calculations show that the lowest-lying scalar, ten- 
sor and axial vector glueballs lie in the mass region of 
1-2.5 GeV. While there is a long history of glueball mass 
calculations in lattice QCD, little is known about the 
glueballs besides their masses. Accurate lattice calcu- 
lations of their size, matrix elements and form factors 
would help considerably in their experimental identifica- 
tion. 

Glueball wave functions and sizes have been studied in 
the past @,0,Hli, 

but much of the early work contains 
uncontrolled systematic errors, most notably from dis- 
cretisation effects. The scalar glueball is particularly sus- 
ceptible to such errors for the Wilson gauge action, due to 
the presence of a critical end point of a line of phase tran- 
sitions in the fundamental-adjoint coupling plane. As 
this critical end-point (which defines the continuum limit 
of a <fi 4 scalar field theory) is neared, the coherence length 
in the scalar channel becomes large, which means that the 
mass gap in this channel becomes small; glueballs in other 
channels seem to be affected very little. Results in which 
the scalar glueball was found to be significantly smaller 
than the tensor were most likely due to contamination of 
the scalar glueball from this non-QCD critical point [7|. 
On the other hand, the calculations using operator over- 
laps obtained from variational optimization for improved 
lattice gauge action, which are designed to avoid spurious 
critical point, show that the scalar and tensor glueballs 
were of typical hadronic dimensions [HQ. A straightfor- 
ward procedure to address the controversy over glueball 
size is to measure the glueball wave function, much in the 
same way as the meson and baryon wave functions were 
measured 10 . 
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In this paper, we study the low-lying scalar and ten- 
sor glueballs and their wave functions with renormalizcd 
tadpole improved Symanzik gauge action [Til ] . Our tech- 
niques for calculating the glueball wave functions from 
Wilson loop operators are outlined in Sec. [Til We present 
and discuss our results at zero temperature in Sec. IIIII 
We extend our method to examine the wave functions 
and masses at finite temperature in this section. Here we 
give an explicit interpretation of deconfincment in terms 
of the power-law behaviour of the correlation function. 
Section IIVI is devoted to the summary and concluding 
remarks. 



II. WAVE FUNCTIONS OF GLUEBALLS 

In contrast with the techniques used in Ref. Q, we 
measure our lattice operators from spatially connected 
Wilson loops. Glueballs are colour singlet states and one 
should be able to construct them with closed-loop paths 
which are gauge invariant. The choice of such loops elim- 
inate the need for gauge fixing 1 . Although the calcula- 
tions in Ref. Q have produced some interesting results, 
the approach suffers from a basic problem: the observ- 
ables are calculated from a lattice version of the 2-glue 
operator, which risks a mixture of glueball states with 
flux states 2 . 

In this study we take a more direct approach to the 
problem. We measure the observables in a three-step 



1 It should be noted that gauge-invariant Wilson loops have a 4 de- 
pendence compared to a 2 dependence of two-link operator used 
in Ref. Q. The lower dimension operator yields a linear depen- 
dence in the correlation function as opposed to a 5 dependence for 
Wilson loops. This improves the glueball signal as a is reduced. 

2 The link-link operator used in Ref. @] sums up a large number 
of loops; some of these loops have a zero winding number and 
project on glueballs - others have non-zero winding number and 
project on flux states also called torelons. 
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procedure. First, we calculate the lattice operator 



$(f,t) =$^[0(£,t) 



(1) 



where is the plaquette operator and $ measures the 
two plaquette or two-loop component of the glueball wave 
function. The r dependence will be reflected in the length 
of links required to close the loops. From a suitable lin- 
ear combinations of rotation, parity inversions and real 
or imaginary parts of the operators involved in <J>, one 
can construct glue ball operators with desired quantum 
numbers @, fl2lll3l| . Since we want to explore the nature 
of wave functions, we focus only on the low-lying "sym- 
metric" and "antisymmetric" scalar channels (which are 
the cosine and sine, respectively of the Wilson loop in 
question) and tensor glueball states. 

The wave function and mass are obtained from the 
correlation function: 



C(f,t) = (^(^$(0,0)), 



(2) 



where one needs to subtract the vacuum contribution 
from the correlator for ++ state. The source can be 
held fixed while the sink takes on the r dependence. This 
proves to be helpful in maintaining a good signal. The 
disentangling of the glueball and torelon is usually taken 
care of automatically by the choice of Wilson or Polyakov 
loops. 

To increase the overlap with the lowest state and re- 
duce the contamination from higher states, we exploit 
the APE link smearing techniques 14j. The procedure 



is implemented by an iterative replacement of the origi- 
nal spatial link variable by a smeared link. This results 
in correlations which reach their asymptotic behaviour at 
small time separations. In addition, the noise from ultra- 
violet fluctuations is reduced. The smearing parameter 
is fixed to 0.7 and ten iterations of the smearing process 
are used. To find the optimum smearing value, n, we 
examine the ratio (at r — and 1) 

C(r,t + 1)/C(r,t), 

which should be maximum for good ground state domi- 
nance. Using lxl loop as template, the best signal is 
obtained with four smearing steps, with lxl and 2x2 
loops being almost indistinguishable. At (3 — 2.0, the 
signal in 1 x 1 showed a slow convergence with n, hence 
2x2 loops were preferred for optimum overlap. A typi- 
cal value which proved to be sufficient for this case was 
n = 2. 

A second pass was made to measure the optimized cor- 
relation matrices 

C t3 {t) = (SCn, 0)) - ($(n))($(r J -)). (3) 

Let ip( k > be the radial wave function of the fc-th eigen- 
state of the transfer matrix, then 



(4) 



The glueball masses and the wave functions are extracted 
from the Monte Carlo average of (t) by diagonalizing 
the correlation matrices C(t) for successive times t: 



C(t) = R(t)D(t)R(t), 



(5) 



where D is a diagonal matrix of the eigenvalues and R 
a rotation matrix whose columns are the eigenvectors of 
C. Each eigenvector of C matches an eigenstate ip( k > (r) 
of the complete transfer matrix. As the wave function 
is largest at the origin, one would first determine the 
glueball mass with the optimal separation, and then fix 
that mass for all r, and extract the wave function for less 
optimal separations. Similar to the case of mesons [l5T |. 
the wave function is expected to decrease exponentially 
with the r at large separations and is therefore fitted with 
the simple form 



(6) 



to determine the effective radius r$. The effective mass 
can be read off directly from the largest eigenvalue cor- 
responding to the lowest energy 



l eff 



= log 



Ao(r = 0,*=l) 



A (r = 0,t = 2) 



(J) 



III. SIMULATION RESULTS AND DISCUSSION 
A. Results at zero temperature 

Most of our Monte Carlo calculations are carried out 
on 16 2 x 16 lattice with periodic boundary conditions (16 2 
is the space-like box and 16 is the extension in Euclidean 
time direction). The gauge configurations are generated 
using the Metropolis algorithm. After the equilibration, 
configurations are stored every 250 sweeps; 3000 stored 
configurations are used in the measurement of glueball 
masses. Measurements made on the stored configura- 
tions are binned into 10 blocks with each block contain- 
ing an average of 300 measurements. The mean and the 
final errors are obtained using single-elimination jack- 
knife method with each bin regarded as an independent 
data point. Three sets of measurements were taken at 
f3 = 2.0, 2.25 and 2.5. Some finite-size consistency checks 
are done at [3 = 2.25 on an 20 2 x 20 lattice. 

The glueball correlation function for the ++ channel 
against t at ft = 2.25 is shown in Fig. [1] It can be 
seen that the expected behaviour of correlation function 
is attained virtually straight away. The absolute errors 
in the correlation functions are expected to be indepen- 
dent of t for large t. Our errors are consistent with this 
expectation. 

Effective mass plot for (3 = 2.5 simulation is presented 
in Fig. [5] For ++ and channels each it was possi- 
ble to find a fit region i mm — imax in which convincing 
plateaus were observed. The effective masses are found 
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lines. The scalar wave function shows the expected be- 
haviour for all the (3 values analysed here. As for the 



FIG. 1: Correlation function for the ++ channel against t. 
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stable using different values of t in Eq. 0, which sug- 
gests that the glueball ground state is correctly projected. 
At (3 = 2.5, we noticed considerable fluctuations in the 
tensor mass at large t. An acceptable fit was only possi- 
ble for t = [2 — 5]. To ensure the validity of our results, 
we compared them to those obtained using 



l eff 



log 



A (t-1)-A (t) 
Ao(t)-Ao(t + l) 



(8) 



It was found that the evaluations of Eqs. ([7]) and 
yielded very consistent within statistical errors. 




FIG. 2: Effective mass plot for scalar and tensor glueball 
states for j3 = 2.5. The dashed horizontal lines indicate the 
plateau values. 

The Wave functions are extracted at time-separations 
t = 1 and 2. We found a little variation (less than two 
percent) in the eigenvectors of C(t) with t which suggests 
that there is no mixing with states of distinct masses. 
Typical plots of the wave functions, normalised to unity 
at the origin, for the symmetric and antisymmetric scalar 
glueballs, at f3 = 2.0, 2.25 and 2.5 are shown in Figs. [3] 
and 2] respectively. For guiding the eyes the Monte Carlo 
points of the same /3-value are connected with straight 



FIG. 3: Scalar ++ glueball wave functions measured on a 
16 3 lattice for various values of (3. 

antisymmetric channel we notice the presence of nega- 
tive contributions in the glueball wave function for r > 6 
at /3 = 2.5. However these contributions do not persist 
when the lattice size is increased from L = 1 6 to 20 (Fig. 
[5]). This would mean that these effects are unphysical 
and can be described as a finite volume artifact. 

For this reason we extract the effective radius of the 
antisymmetric state from the wave function obtained at 
larger volume 3 . The symmetric scalar glueball wave func- 
tion, on the other hand, barely changes sign. 
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FIG. 4: Scalar glueball wave function on a 16 lattice 
at (3= 2.0,2.25 and 2.5. 



3 The results for 16 3 lattice in Fig. [4] are shown only as an illus- 
tration. Comparison of the data for the effective mass on two 
lattice sizes reveals that none of our states could be interpreted 
as a torelon pairs, since no mass reduction of sufficient magnitude 
was found as the lattice volume was reduced. 
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FIG. 5: Scalar glueball wave function on a 20 3 lattice 
at (3 = 2.25 (open circles) and /3 = 2.5 (open triangles). 




FIG. 6: Tensor glueball wave function on a 16 3 lattice at 
(3 = 2.0,2.25 and 2.5. 



Fig. [6] shows the wave function for the tensor glueball, 
at (3 = 2.0, 2.25 and 2.5. The tensor wave function re- 
mains positive and shows the expected flatness. It can 
be seen that tensor glueball is much more extended than 
the scalar as one moves towards higher (3 values. This 
would imply that tensor is therefore more sensitive to 
the finite-size effects, which is very visible in the distor- 
tion of the wave function for large r at f3 = 2.5. Naively 
we would expect that the spatial size at which we begin 
to encounter large finite size effects to be related to the 
size of the glueball. 

The expected finite-size scaling behaviour of the mass 
gap near the continuum critical point in this model is 
not known; but Weigel and Janke [IB] have performed 
a Monte Carlo simulation for an 0(2) spin model in 
three dimensions which should lie in the same univer- 
sality class, obtaining 



M fa 1.3218/L 



(9) 



FIG. 7: The logarithmic plot of the scalar and tensor glueball 
wave functions at j3 = 2.25 and 2.5. The effective radii can 
be obtained from the inverse slopes of the curves. 



tions on a 20 2 x 20 lattice at (3 = 2.25 and 2.5. The mass 
and size of ++ channel are almost unchanged as the lat- 
tice size increases from 16 to 20. We also find that our 
estimates for the tensor state are consistent with no fi- 
nite volume dependence at (3 = 2.25. However, the tensor 
mass was found to increase by about 4% and the effective 
radius by about 7% from 16 to 20 lattices at (3 = 2.5. We 
do not have enough data extrapolate mass and the radius 
to the infinite volume limit or to check whether the dif- 
ference is due to statistical errors or whether there is an 
incomplete convergence. Given that no mass reductions 
of sufficient magnitudes were found as the lattice volume 
is changed, none of our states could be interpreted as a 
torelon pair. 

In order to get some quantitative information on the ef- 
fective radius, the glueball wave functions are fitted in the 
range 3 < r < 8 by the form ([B]). This form fits the data 
rather well for the scalar glueball with the best-fit esti- 
mates obtained with a x 2 /Ndf of 0.92 - 0.67. Due to dis- 
tortion 4 in the tensor wave function at small r at /3 = 2.5, 
a meaningful fit was possible only in the range 6 < r < 8. 
The effective radius obtained was confirmed by examin- 
ing the plateau in the ratio \og\ip (r) /ip(r + 1)]. Note, 
that our lograthmically plotted wave functions (Fig. [7]) 
are merely illustrations. 

To summarize: in the weak coupling region a spectrum 
of massive ++ , 0~~ and 2 ++ glueballs is indicated with 

m(0~~ ) < m(0 ++ ) < m(2 ++ ). 

Since there is a good signal for wave function persisting 
long enough to demonstrate convergence to the asymp- 
totic value, it seems to be reasonable to estimate mass 



for the magnetic gap. In order to ascertain the finite-size 
effect on our measurements, we performed extra simula- 



4 Because of the distortion and impossible complete elimination of 
all the excited-states, especially near r ~ 0, it follows that Eq. 
jE\ holds only in the limited interval, which does not include the 
vicinity of r ~ 0. 



ratios with our present method. The estimates of masses 
and ro, in lattice units, at various (3 values are shown in 

Tables Hn and EH 



TABLE I: Masses of scalar glueballs in lattice units for two 
spatial extensions, L — 16 and 20. 



Mass 



0++ 



/3/L 16 


20 


16 


20 


2.0 0.803(6) 




0.441(4) 




2.25 0.523(3) 


0.527(3) 


0.266(3) 


0.261(4) 


2.5 0.364(3) 


0.369(2) 




0.182(2) 


TABLE II: Sizes of scalar glueballs in lattice units for two 


spatial extensions, L - 


= 16 and 20. 








Size 






0++ 




0" 




/3/L 16 


20 


16 


20 


2.0 1.4(2) 




1.0(1) 




2.25 2.75(4) 


2.8(3) 


2.14(4) 


2.2(2) 


2.5 5.1(1.0) 


5.2(9) 




5.0(1.0) 



TABLE III: Mass and size of tensor glueballs in lattice units 
for two spatial extensions, L — 16 and 20. 





Mass 






Size 


/3/L 


16 


20 


16 


20 


2.0 


1.2(1) 




5.0(7) 




2.25 


0.82(2) 


0.81(6) 


9.7(1.7) 


9.8(1.4) 


2.5 


0.544(2) 


0.58(1) 


10.1(2.6) 


10.2(2.4) 



Our results for lattice masses and mass ratios are gen- 
erally, within statistical errorSjjn agreement with the ex- 
isting Euclidean estimates 

Bill E3, if Perhaps a little 
high in places. Qualitatively our results, at zero temper- 
ature, are in agreement with the scenario of spectrum of 
massive magnetic monopoles. 



TABLE IV: Glueball sizes in the units of string tension. 



a 


K(= a 2 a) 






r —y/a 




2.0 


0.0508(5) 


0.0856 


0.31(14) 


0.24(9) 


1.13(18) 


2.25 


0.0221(3) 


0.0481 


0.40(17) 


0.32(16) 


1.14(22) 


2.5 


0.0119 


0.0272 


0.56(21) 


0.50(19) 


1.11(25) 



To extrapolate our effective radii to the continuum 
limit, we take the dimensionlcss products of sizes so that 
the scale, a, in which they are expressed cancels. We 
choose to take products of the effective radii, rop/a, to 



a^fa since the string tension is our most accurately cal- 
culated quantity. As in the (3+l)D confining theories, 
we expect that dimensionless product of physical quanti- 
ties, such as ro^v^j wm approach their continuum limit 
with correction of O(a^i), where a e // is the effective lat- 
tice spacing in "physical units" when the mass gap has 
been renormalized to a constant The string tension, 
K{= a 2 er), is obtained by using the Wilson loop aver- 
ages and fitting the on-axis data with V(r). In Fig. [8] 
we show the product r^^^fa plotted against a e ff- Since 
the products are plotted against a e f f , the continuum ex- 
trapolations are simple straight lines. We notice that the 
product ro t 2*/a varies only slightly over the fitting range. 
The non-zero lattice spacing values of the product are 
within 0.04 - 0.29 and 0.01 - 0.02 standard deviations of 
the extrapolated zero lattice spacing results for the scalar 
and tensor glueballs respectively. The striking feature of 
this plot is the little variation of the product with a e f / . 
This will make for very accurate and reliable continuum 
extrapolations. Linear extrapolations to the continuum 
limit yield values of 0.60 ± 0.05 and 1.12 ± 0.03, in the 
units of string tension, for the scalar and tensor states, 
respectively. In contrast to the tensor, the scalar glueball 
size shows significant finite-spacing errors. By setting the 
string tension to 420 MeV, we obtain the physical radii 
of 0.28(7) and 0.52(5) fm, for the scalar and tensor glue- 
balls, respectively. Our results show the size of the tensor 
glueball roughly two times as large as the scalar glueball. 
These estimates agree with the rough estimates of glue- 
ball sizes obtained at various temperatures in Ref. 0]. 
This is an improvement over the estimates obtained in 
Ref. 0] where the predicted radius for the tensor glue- 
ball (~ 0.8 fm) was found four times larger than scalar 
glueball radius. 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 



FIG. 8: Glueball radii in the units of string tension as a 
function of the effective spacing, a e //. Extrapolations to the 
continuum limit are shown as dashed lines. 
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B. Finite temperature results 

To check the consistency of our method, we performed 
a study on an asymmetric lattice: 16 2 x 4 at f3 = 2.25. 
The procedure itself is a straight forward extension of the 
procedure adopted in the previous subsection. We do not 
plan to study the high temperature aspects of this model 
here but focus on the behaviour of the glueball mass and 
wave function in the deconfined region. 

The physical temperature T = l/(aN t ), is given via 
the lattice parameters as follows: 



1 



N t VK 



(10) 



For completeness, we give a temperature estimate of 
1.125 in the units of string tension. By setting the string 
tension to 420 MeV, we estimate a physical temperature 
of T ~ 1.25T C , where the T c ~ 360 MeV at pseudo- 
critical coupling /3 C = 1.87(2) 18]. One expects [H 
that the high temperature phase has a massless photon 
and the linear potential is replaced by the two dimen- 
sional logarithmic Coulomb potential. This logarithmic 
behaviour is equivalent to a power-law dependence of the 
Wilson loop correlation function, 



C(r) = (pt(r)P(0)) 



-V(T) 



(11) 



with an exponent which decreases as T increases. Fur- 
thermore, since the high-temperature phase of the gauge 
theory corresponds to the ordered phase of the spin sys- 
tem, the predicted power-law behaviour of the correla- 
tion function is just like that of a two-dimensional U(l)- 
invariant spin system - a 2-D XY model. 

Fig. [9] shows a plot of correlation functions versus sep- 
aration. The straight line indicates the fit to the form 
(jlip . The finite temperature phase transition is visible 




log(r) 



FIG. 9: The logarithmic plot of the correlation function at 
[3 — 2.25. The straight line indicates power-law behaviour. 

in the change of the correlation function from exponen- 
tial to power-law behaviour. Thus it becomes evident 



that T > T c in our simulation. It can be seen that form 
(fTl) fits the data rather well. Nonetheless, our Monte 
Carlo simulations were unable to confirm that the expo- 
nent is mo ving towards the value of 0.25 (that of the 2-D 
XY model [20() predicted for the continuum theory. Our 
estimated value for the exponent is four times larger than 
the predicted value. This indicates that our P value of 
2.25 is not large enough to gives us reason to hope that 
we are approaching continuum physics. An interesting 
feature to explore in this context is whether the cou- 
pling to the matter fields in the leading order (/? — ► oo) 
calculations will move the critical exponent towards the 
predicted value. 
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FIG. 10: Scalar and tensor glueball wave functions measured 
on 16 2 x 4 lattice at /? = 2.25. 

In the deconfinement phase above the critical temper- 
ature, glueballs are no more elementary excitations. At 
high temperatures we have a plasma that behaves in bulk 
roughly like a free gas of quarks and gluons, thus forming 
a new phase, i.e., the quark-gluon plasma (QGP) phase. 
Above the critical temperature T c , properties such as 
confinement and chiral symmetry breaking disappear. A 
detailed understanding of thermal glueballs gained over 
the last decade can be found in [§, [2l|, [2^, [H, Q, [25[ and 
the references therein. As a result, quarks and gluons 
are liberated and tremendous changes are expected in the 
mass spectrum. Fig. [TU]shows the scalar and tensor wave 
functions obtained through the same analysis as in Figs. 
[3] and [H Our results indicate that no significant changes 
occur in the scalar and tensor wave functions. Glueball 
masses appear almost unaffected. By comparing the re- 
sults at T — and T = 360, we observe an effective mass 
reduction, {am G (T - 0) - am G (T ~ 360)), of about 4%, 
with statistical uncertainties typically on less than a per- 
cent level, for ++ and 2 ++ glueball modes. This appears 
to be a very small change since we expect a rather contin- 
uous mass reduction of glueballs in the deconfined phase. 
This might be due to the fact that for zero momentum the 
power-law behaviour of the correlation function leads at 
short distances to the spin-wave results, which prevents 
us from seeing the massless excitations. 
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The non-vanishing effective masses would suggest the 
presence of glueball modes above T c . Other work on fi- 
nite temperature SU(3) 0, Hl| has also confirmed the 
survival of correlations above T c in the scalar and tensor 
colour-singlet modes. However, these studies have shown 
that thermal mass changes rather continuously across the 
critical temperature. The existence of the effective mass 
gives rise to the possibility that some of the nonpertur- 
bative effects survive in the deconfined phase, and the 
colour-singlet modes exist as metastable states above T c . 
The metastable states in the ordered phase (large j3) ap- 
pear to be caused by the unusually large separation of a 
vortex pair, which may take many sweeps to recombinc. 
Near the transition the number of vortices increase, and 
some of them begin to unbind. This eventually drives 
the system into a disordered phase as one moves to the 
region T < T c . 

Whether bound or metastable modes, the glueballs can 
decay into two or more gluons thus acquiring finite width 
which is expected to become less negligible in the decon- 
fined phase. Thus it becomes important to take into ac- 
count the effect of width in best-fit analysis. This might 
also explain a very modest reduction of our masses at 
T > T c . However, from this study, it is not possible 
to determine whether such colour singlet modes really 
survive above T c as metastable states. An extensive sys- 
tematic analysis, of unquenched improved lattice QCD 
at finite temperature, along these lines is under way [261 ] - 

IV. SUMMARY AND CONCLUSION 

We have studied wave functions and sizes of scalar and 
tensor glueballs using improved 3-dimensional U(l) lat- 
tice model. In this preliminary study we take a more 
direct approach to the problem; instead of fixing a gauge 
or a path for the gluons, we measure the correlation 
functions from our lattice operators from spatially con- 
nected Wilson loops which, being the expectation values 
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of closed-loop paths, are gauge invariant. This approach 
has the advantage that the disentangling of the glueball 
and torelon is usually taken care of automatically by the 
choice of Wilson or Polyakov loops. We observed that 
the size of tensor glueball is roughly two times larger 
than the size of the scalar glueball. We believe that our 
estimates are more reliable than the results obtained in 
Ref. where the size of the tensor glueball was found 
to be ~ 0.8 fm, four times as large as its scalar counter- 
part. The predicted zero lattice spacing results are not 
actually found by extrapolation to zero lattice spacing, 
but are obtained instead from calculations at (3 of 2.2 of 
glueball size, with no accurate representation of the effect 
of the absence of extrapolation. Also the results were of 
limited interest because of their manifest dependence on 
the gauge chosen and the problem of disentangling of the 
glueballs and torelons. 

Finally, for completeness, we extended our method to 
measure the wave function and mass for a finite temper- 
ature deconfinement phase. For the lowest ++ and 2 ++ 
glueballs, no significant mass reduction was observed in 
the deconfined phase, while the wave functions remain 
almost unchanged. The existence of the effective mass 
indicates that colour-singlet modes may survive in the 
deconfined phase as metastable states. In such a case 
glueball decay and decay width, as spectral component, 
in the deconfinement phase are the most feasible candi- 
dates for a more reliable analysis for the future studies. 
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